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ABSTRACT 


The  theory  of  gain  saturation  of  a  two -level  system  is  briefly  reviewed  and  then 
generalized  to  include  an  arbitrary  number  of  upper  and  lower  levels;  the  upper 
levels  are  coupled  to  each  other,  as  are  the  lower  levels.  The  laser  radiation 
is  assumed  to  induce  transitions  between  one  upper  and  one  lower  level.  In  the 
limit  of  tight  coupling  among  the  upper  and  lower  levels  by  themselves,  the 
equations  reduce  to  those  of  a  two-level  system,  with  effective  relaxation  rates 
that  are  weighted  sums  of  the  relaxation  rates  of  the  multilevel  system.  The 
relaxation  rates  occurring  in  a  C07  laser  system  are  comparatively  low  so  that 
one  would  expect  spatial  diffusion  to  play  an  important  role  in  determining  sat¬ 
uration.  A  theory  of  diffusion  is  carried  through  and  it  is  shown  that  diffusion 
effects  can  indeed  be  important  for  optical  beam  diameters  of  a  few  millimeters . 
Finally,  experiments  on  a  sealed -off  CC^  laser  oscillator  and  amplifier  system 
are  reported.  Amplifier  gain  is  measured  for  four  different  beam  radii.  The 

"equivalent"  saturation  parameter  derived  from  the  measurements  decreased 

2 

monotonically  from  97  to  25  W/em  as  the  average  input  beam  radius  increased 
from  0.  9  to  2.5  mm  in  the  9 -mm  radius  discharge  tube  of  the  amplifier. 


Accepted  for  the  Air  Force 

Franklin  C.  Hudson 
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PREFACE 


The  bulk  of  the  work  presented  in  this  report  has  appeared  as  a  Master’s  Thesis 
by  C.  P.  Christensen,  submitted  to  the  Department  of  Electrical  Engineering  at 
Massachusetts  Institute  of  Technology  in  September  1968,  entitled  "Diffusion, 
Relaxation  and  Gain  Saturation  in  the  CO^  Laser."  Although  based  on  the  thesis, 
this  report  differs  in  several  ways.  When  the  thesis  was  composed,  we  were  un¬ 
certain  to  what  extent  the  other  rotational-vibrational  levels  (e.g.,  the  J-levels 
of  00° 2)  could  contribute  to  the  "effective  number  of  levels."  Further,  the  num¬ 
bers  used  in  the  thesis  to  estimate  the  importance  of  diffusion  were  very  rough. 
More  detailed  work  by  H.  A.  llaus  gave  more  accurate  magnitudes.*  Sections  I 
through  III  were  accomplished  with  close  cooperation  of  Christensen  and  Haus; 
Sec.  IV  is  the  work  of  H.  A.  llaus,  and  Sec.  V  that  of  C.  P.  Christensen.  The  ex¬ 
periments  of  Sec.  VI  were  conducted  at  Lincoln  Laboratory  by  C.  P.  Christensen 
and  C.  Freed.  An  abbreviated  version  of  this  work  centered  around  the  exper¬ 
imental  results,  under  the  title  "Gain  Saturation  and  Diffusion  in  CO^  Lasers"  by 
C.  P.  Christensen,  C.  Freed,  and  H.  A.  Haus,  will  be  published  in  the  June  1969 
issue  of  the  IEEE  Journal  of  Quantum  Electronics. 

H.  A.  Haus 


*  This  work  was  supported  by  Joint  Services  Electronics  Program  (Contract 
DA28-043-AMC  02536(E)). 


EFFECT  OF  DIFFUSION  ON  GAIN  SATURATION 
IN  C02  LASERS 


INTRODUCTION 

In  a  laser  system,  the  gain-saturation  parameter  determines  how  the  gain  decreases  with 

increasing  intensity.  The  CO?  laser  is  a  high-power  laser,  in  part,  because  it  has  a  large  gain- 
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saturation  parameter.  Published  experimental  results  on  the  CO^  laser  reported  rather 
widely  different  saturation  parameters.  The  conditions  in  these  experiments  were  sufficiently 
different  so  that  the  variances  in  the  results  in  themselves  would  not  be  too  surprising.  Yet,  the 
work  reported  here  was  stimulated  in  part  by  the  rather  large  variation  in  the  saturation  param¬ 
eter  found  by  these  workers,  in  the  belief  that  there  existed  some  unsuspected  cause  for  the  dis¬ 
crepancies.  The  desire  to  predict  the  power  emitted  by  a  CO^  laser  at  threshold  was  another 
reason  for  the  investigations  reported  here;  this  prediction  requires  knowledge  of  the  saturation 
parameter.  The  saturation  parameter  predicted  on  the  basis  of  a  two-level  description  of  the 
CO^  laser  system  is  an  order  of  magnitude  too  low  compared  with  the  experimentally  observed 
value.  Hence,  it  is  immediately  apparent  that  a  multilevel  description  is  called  for,  in  view  of 
the  fact  that  the  coupling  among  the  vibrational-rotational  levels  of  the  00°  1  and  10 °0  states  is 
strong.  Further  study  led  to  the  realization  that  diffusion  effects  must  play  an  important  role 
in  determining  the  gain  saturation  in  a  CO^  laser  with  a  beam  radius  of  up  to  a  few  millimeters. 

In  Sec.  I,  we  present  a  brief  review  of  gain  saturation  in  a  two-level  system.  Although  the 
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results  are  known  and  have  been  used  in  the  literature,  the  original  derivation  has  never 
been  published.  In  presenting  it  here,  we  provide  a  useful  background  for  the  discussions  in 
Sec.  II  of  the  multilevel  system  with  a  lasing  transition  occurring  between  two  levels,  still  dis¬ 
regarding  spatial  diffusion.  In  Sec.  Ill,  we  specialize  to  the  case  of  a  multilevel  system  with 
tight  coupling  among  its  upper  and  lower  levels.  A  closed-form  expression  can  be  obtained  for 
the  gain-saturation  parameter  depending  only  upon  the  Boltzmann  equilibrium  densities  and  the 
relatively  low  relaxation  rates  of  the  levels  considered  into  the  levels  not  explicitly  included  in 
the  rate  equations. 

In  Sec.  IV,  we  discuss  diffusion  in  the  mixture  of  gases:  CO^,  N^,  He,  and  A  diffusion 

constant  is  obtained  for  the  pressures  and  temperature  used  in  a  typical  sealed-off  CO^  system. 
We  find  that  a  CO^  molecule  can  diffuse  across  an  optical  beam  of  a  few  millimeters  in  a  time 
comparable  to  the  inverse  relaxation  rate  of  the  lasing  level;  hence,  we  conclude  that  spatial 
diffusion  effects  can  play  an  important  role  in  determining  gain  saturation  in  C O^.  Section  V  is 
an  analysis  of  gain  saturation  in  the  presence  of  spatial  diffusion  for  an  optical-beam  profile  of 
rectangular  shape.  This  analysis  can  be  carried  out  in  closed  form,  and  gives  an  estimate  of 
the  influence  of  diffusion  upon  the  gain-saturation  parameter.  The  experimental  results  are 
presented  in  Sec.  VI. 
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I.  THE  TWO- LEVEL  SYSTEM 


We  first  investigate  gain  saturation  in  the  steady  state  in  a  two-level  system.  A  rate  equa¬ 
tion  approach  is  adequate  because  of  the  steady-state  condition,  and  spontaneous  emission  can 
be  disregarded  because  the  power  levels  to  be  considered  are  high.  We  have  the  rate  equations 


dN  XT 

— -  —  —  y  N  — W 

dt  'u 
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dt 


=  —  y  £  n  —  W 
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(1-2) 


where  N  is  the  population  density  of  the  upper  level,  n  is  that  of  the  lower  level,  the  y's  are  the 
relaxation  rates,  and  W  is  the  rate  of  induced  emission  which  is  related  to  the  spontaneous  tran- 
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I  =  intensity  of  radiation, 

X  =  wavelength  of  radiation, 
v  =  frequency  of  radiation. 


g^) 
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with 


v  -  center  frequency  of  lasing  transition. 


Av  -  Lorentzian  half -width  of  transition; 

gu  and  g^  are  the  degeneracy  factors  of  upper  and  lower  levels,  and  the  R’s  are  the  pumping 
rates. 

The  gain  constant  a  is  obtained  from  the  above  by  noting  that,  by  definition. 
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on  one  hand,  and,  from  energy  conservation. 
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on  the  other.  Hence, 
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Returning  to  the  rate  equations  and  noting  that  d/dt  0  in  the  steady  state,  we  may  solve 
for  [N  —  (gu/gj j)  n]  in  terms  of  I,  obtaining 
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When  the  radiation  is  weak,  W  in  the  denominator  of  Eq.  (1-8)  can  be  set  equal  to  zero,  and  we 
obtain  the  small-signal  gain  c*o. 


\ru  g,  Vf ) 


ao  8  7rt_ 


sp 


With  increasing  intensity,  the  gain  decreases  according  to  the  law 
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is  the  saturation  parameter  at  the  center  of  the  Lorentzian  line. 

Generally,  1  /y ^  «  l/*yu  and,  in  this  case,  Eq.  (1-11)  can  be  rewritten,  by  returning  to  the 
definition  of  W  in  Eq.  (1-3), 

7U  =  W(I  ~  Is)  •  (1-12) 

Thus,  saturation  occurs  according  to  Eq.  (1-12)  when  the  rate  of  induced  transitions  becomes 
comparable  to  the  relaxation  rate  of  the  upper  level.  The  induced  transitions  depopulate  the  up¬ 
per  level  at  a  rate  comparable  to  the  rate  at  which  molecules  enter  or  leave  the  upper  level  by 
relaxation  processes  in  the  lasing  medium.  Thus,  when  considering  a  system  of  given  small- 
signal  gain,  we  can  increase  the  saturation  parameter  (and  hence,  presumably,  the  output  power 
of  a  laser)  by  increasing  the  relaxation  rate.  At  first  glance,  this  seems  somewhat  paradoxical; 
indeed,  it  seems  that  a  system  is  made  harder  to  saturate,  and  thus  capable  of  delivering  a  higher 
power  density,  when  we  increase  the  rate  at  which  particles  leave  the  excited  level  (and  cease 
to  interact  with  the  lasing  field).  The  paradox  disappears  as  soon  as  we  note  that  by  raising  the 
rate  of  relaxation  of  the  upper  level  at  a  given  small-signal  gain,  we  must  raise  the  pumping 
rate;  in  other  words,  we  must  raise  the  rate  at  which  particles  enter  the  upper  level. 

One  means  of  effectively  raising  the  rate  at  which  particles  start  or  cease  to  interact  with 
the  lasing  field  is  by  diffusing  particles  into  and  out  of  the  laser  beam.  Hence,  diffusion  can  be 
expected  to  increase  the  saturation  intensity.  The  increase  would  become  noticeable  when  the 
rate  at  which  particles  leave  the  beam  via  diffusion  becomes  comparable  to  the  rate  of  relaxation 
of  the  upper  level  y  .  We  shall  return  to  this  point  after  studying  the  multilevel  system. 
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n.  THE  MULTILEVEL  SYSTEM 

We  consider  a  system  consisting  of  a  set  of  (tightly)  coupled  upper  laser  levels  and  a  set  of 
(tightly)  coupled  lower  levels.  The  weaker  coupling  to  all  other  levels  is  represented  by  a  phe¬ 
nomenological  relaxation  rate  to  each  of  the  upper  and  lower  levels.  Lasing  action  is  assumed 
to  occur  only  between  one  upper  and  one  lower  level,  denoted  by  the  subscript  1.  The  lasing 
levels  are  assumed,  at  first,  to  be  homogeneously  broadened.  Spontaneous  transitions  will  be 
neglected. 

In  the  C system,  the  set  of  upper  levels  is  assumed  to  represent  the  tightly  coupled  set 
of  upper  vibrational-rotational  levels  (00 °n,  J),  the  lasing  level  being  the  (00°  1,  J^)  level.  The 
rate  equations  for  this  system  are  of  the  form 
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Here,  and  n^  are  the  population  densities  in  the  upper  and  lower  lasing  levels,  respectively. 
The  remaining  N^'s  and  n^’s  are  the  population  densities  of  the  upper  and  lower  vibrational- 
rotational  levels,  respectively.  R^1  and  R^  are  the  pumping  rates,  and  W  is  the  induced  transi¬ 
tion  rate  defined  in  Eq.  (1-3).  The  T.V  *s  and  *s  are  the  relaxation  rates  among  the  upper  and 

u  ^  J  tli 

lower  levels,  respectively,  y,  is  the  rate  of  relaxation  of  the  k  upper  level  into  all  other  levels 
not  explicitly  included;  y^  is  the  rate  of  relaxation  of  the  kl  1  lower  level  into  all  levels  not  ex¬ 
plicitly  included.  The  remaining  quantities  were  defined  in  Sec.  I. 

In  the  steady  state,  the  time  derivatives  in  Eq.  ( 1 1  - 1 )  are  equal  to  zero.  The  equations  can 
be  solved  for  and  n^  in  terms  of  W(I)  by  means  of  Cramer’s  rule  and,  in  this  way,  we  obtain 
an  expression  for  gain  saturation  as  a  function  of  intensity.  First,  define  the  effective  relaxa¬ 
tion  rate  for  the  k^1  level 
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Defining  the  determinants 
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and  similar  determinants  for  the  quantities  referring  to  the  lower  levels,  we  obtain  for  the  pop¬ 
ulation  density  in  the  upper  lasing  level 
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A  similar  expression  is  obtained  for  the  lower  level. 
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From  Eqs.  (II-6)  and  (II-7),  we  obtain  the  population  inversion.  From  the  fact  that  the  power  added 
to  the  laser  beam  per  unit  volume  is  hr  [N^  —  (gu/g^)  n^]  W,  we  compute  the  gain  constant  a 
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Several  features  are  worth  noting.  The  dependence  of  a  upon  intensity  I  is  the  same  as  that  of 
a  two-level  system  (the  quantity  in  brackets  is  negative).  Further,  for  any  particular  set  of  dis¬ 
charge  conditions,  all  parameters  are  constants  in  terms  of  which  the  saturation  parameter  (i.e., 
the  intensity  at  which  the  gain  constant  drops  to  half  its  value)  can  be  evaluated 
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In  terms  of  the  saturation  parameter,  we  may  write  the  dependence  of  ca  upon  I: 
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If  the  system  is  inhomogeneously  broadened 


9,10 


the  gain  constant  is  produced  by  interaction  with 


molecules  of  varying  center  frequency  v Q.  The  probability  distribution  P(^Q)  over  v  for  a  gas 
in  thermal  equilibrium  is  Gaussian 
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exp 
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where  v is  the  line  center  frequency  of  a  molecule  with  zero  velocity.  The  gain  constant,  for 
the  case  when  the  Gaussian  linewidth  Ar^  is  much  larger  than  Ar,  is 
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If  we  study  the  decrease  of  n  with  intensity,  at  small  intensity,  we  find  a  linear  decrease  with  I 


(11-11) 


where  I  is  given  by  the  same  expression  as  for  the  homogeneously  broadened  line.  The  gain  de¬ 
creases  less  rapidly  with  intensity  because,  in  effect,  a  larger  population  density  is  participating 


6 


in  the  interaction  than  would  have  been  if  all  molecules  had  the  same  center  frequency  and  had 
a  density  such  that  the  same  small-signal  gain  would  have  been  achieved. 

When  the  inhomogeneous  broadening  is  not  as  large  as  assumed  to  arrive  at  Eq.  (11-10),  the 
integration  cannot  be  carried  out  in  closed  form.  Still,  the  gain  decreases  with  increasing  I; 
initially,  at  small  I,  the  dependence  is  as  in  Eq.  (11-11)  with  the  factor  in  front  of  l/lg  varying  be¬ 
tween  1  and  depending  upon  the  degree  of  inhomogeneous  broadening.  We  may  take  advantage 
of  this  fact,  and  define  a  saturation  parameter  I  by  fitting  the  initial  dependence  of  a  upon  I  to 
the  straight-line  dependence 

o  =  »0  r_f  (n-12) 


where  |  Ig^  Iq^  I  ,  depending  upon  the  degree  of  inhomogeneous  broadening.  The  important 
fact  to  note  is  that  I  is  a  function  of  the  relaxation  rates,  linewidth,  and  degree  of  inhomogene¬ 
ous  broadening,  but  is  not  a  function  of  the  beam  geometry  (beam  diameter)  if  the  rate  equation 

description  Eq.  ( II  —  1 )  is  indeed  the  proper  model  of  the  system. 

u  u  ft  f 

In  Sec.  Ill,  we  shall  make  the  assumption  that  y^  «  and  y^  «  I\^,  and  obtain  an  esti¬ 
mate  for  the  saturation  parameter  I  . 


m.  EQUIVALENT  RELAXATION  RATE 


When  the  coupling  among  the  upper  levels  is  tight  (y «  r.Y  )  and  the  same  holds  for  the 

J  u  f 

lower  levels,  we  can  expand  the  determinants  in  Eq.  (11-11)  in  powers  of  y^  and  y^  ,  and  express 

the  result  in  terms  of  equivalent  relaxation  rates  y  U  and  y  *  . 

u  i  ec*  eC*  u  f 

The  determinants  A  and  A  are  zero  to  zeroth  order  in  y^  and  y^  ,  respectively.  The  first- 
u  a  u  n 

order  term  in  y^  and  y^  of  A  and  A  can  be  found  in  terms  of  the  cofactors  of  the  y's  in  A.  In¬ 
deed,  to  first  order  in  the  y's,  we  may  set  all  but  one  of  the  y's  equal  to  zero  in  A  of  Eq.  ( II -3 ) 
and  then  sum  over  all  the  terms  obtained  in  this  way.  Further  simplification  is  obtained  through 
the  principle  of  detailed  balance. 

Defining, 


y.V,  =  N„er® 

jk  k  jk 


where  N^  represents  the  equilibrium  density,  we  have  from  the  principle 


u  u 
yjk  =  Ykj 


(in-i) 


(III-2) 


and  a  similar  expression  for  the  lower  levels. 

Now  consider  the  cofactor  of  the  11  term  in  the  determinant  Eq.  ( II  —  3 ),  with  all  yu,s  set  equal 
to  zero.  We  find 
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Now  consider  the  cofactor  of  the  22 -element 
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(III  —4) 


Adding  all  columns  to  the  first  column  and  all  rows  to  the  first  row  in  the  determinant 
Eq.  (III-4),  we  find  that 
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(ITI-5) 


Similarly,  we  can  prove  that 


C . .  C,  , 
JJ  _  kk 
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( III-6 ) 


To  lowest  order  in  the  y.  's, 
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Au  =  C 


11 


and 
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M 

2  y  UN  e 
i=i  1  1 


(III-7) 


(III  —8 ) 


i.  i. 

A  comparable  expression  is  found  for  A  /A  .  Introducing  these  into  Eq.  ( II  —  9) ,  we  have  for 
the  saturation  parameter  of  a  homogeneously  broadened  medium 


87rt  hr 
sp 
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One  interesting  observation  may  be  made  immediately:  When  all  y  ’s  and  all  N.  ’s  are 

equal,  then  l/yU  of  a  two-level  system  is  replaced  by  l/My  ,  where  M  is  the  number  of  levels. 

u  l 

A  corresponding  relation  holds  for  the  lower  level.  The  rates  y  and  y  are  increased  by 
a  factor  equal  to  the  number  of  levels,  and  the  saturation  parameter  is  raised  accordingly. 

More  generally,  it  may  be  stated  that  the  gain  saturation  of  the  system  behaves  similar  to 
that  of  a  two-level  system  with  equivalent  relaxation  rates 
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(III-10) 
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We  may  estimate  these  rates  for  the  C02  system  by  assuming  that  all  yU,s  (and  the  yl ’s) 


are  equal  to  each  other.  The  populations  of  the  J 
(2 J  +  1 )  exp 


.th 


rotational  levels  are  proportional  to 


h  c 

-BJ(J  +  1)  “,] 


where  the  parameter  B  is  characteristic  of  the  vibrational  state  to  which  this  vibrational- 
rotational  level  belongs.  The  populations  of  the  vibrational  states  are  in  the  ratio's  of 
exp  [  —  (E/kT )],  where  E  is  the  energy  of  the  vibrational  excitation.  Table  I  shows  the  parameters 
obtained  from  the  literature,  and  from  which  the  effective  enhancement  of  the  relaxation  rate 


TABLE  1 

Vibrational  State 

B 

(cm  ’) 

Energy 
(cm  ') 

Computed  No. 
of  J  Levels 

Relative  Total 
Population  Density 

00°  1 

0.3866 

2349 

~  14 

i 

00°2 

0.3837 

4679 

~  14 

0.033 

00°  3 

0. 3808 

6976 

~14 

0.011 

00°4 

0. 3779 

9256 

~14 

0.007 

can  be  obtained.  We  assume,  in  the  absence  of  more  precise  information,  that  all  yll,s  are 
equal  to  each  other.  The  sum  over  the  rotational  levels  within  one  single  vibrational  state  is 
approximated  by  an  integral 

Yj  (2J  +  1)  exp[-BJ(J  +  i) 

J  even 
or  odd 


foo 

dJ(2J  +  l)exp[  —  BJ(J  +  1,^]=  ^  . 

o 

If  the  lasing  level  is  assumed  to  have  the  largest  population, 

=  0 


jj  (2 J  +  1)  exp  [  — BJ (J  +  i)  ££] 


J=J, 


we  have,  solving  for  and  taking  the  value  B  from  Table  I,  at  T  =  300 °K: 

Ji-Jz 


kT 

r  2Bhc 


-1-18 


(III-12) 
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The  population  density  is  proportional  to  (2J1  +  1)  expf—B^J^J^  +  1)  (hc/kT)].  Summed 
over  all  rotational  states,  the  population  density  in  each  of  the  vibrational  states  is  given  by 
Eq.  ( III-12 ).  Hence,  the  effective  relaxation  time  of  the  upper  level  is 

2  expf-lh^/kTJl/B. 

Yeff  =  yU  Ni6/Nie )"yU  2*kT  (2J1  +  T)  exp[-B1J1(J1  +  1)  (hc/kT)]  '  (IH-13) 

A  similar  expression  is  obtained  for  noting  only  that  the  J  value  of  the  lower  level  in  the 

P-branch  is  higher  by  one  than  that  of  the  upper  level.  Introducing  the  values  of  Table  I,  we 
find  2  ISLe/N^  =  14.1  at  T  =  300 °K.  If  the  "temperature"  of  the  vibrational  distribution  is  larger 
than  that  of  the  rotational  distribution,  then  the  exponentials  exp  [— (hr  ./kT)]  are  larger,  leading 
to  larger  entries  in  the  last  column  of  Table  I. 

IV.  DIFFUSION  IN  A  MULTI-SPECIES  SYSTEM 

The  process  of  gain  saturation  in  a  sealed-off  CO^  laser  under  the  influence  of  a  laser  field 
with  a  small  optical  beam  radius,  such  that  diffusion  plays  a  role,  proceeds  roughly  as  follows. 

The  laser  field  depletes  the  population  in  the  upper  vibrational-rotational  levels  and  causes 
an  excess  population  in  the  lower  levels.  Therefore,  excited  CO^  molecules  from  the  outside  of 
the  beam  region  diffuse  into  the  "spatial  hole"  formed  in  the  population  density  of  the  upper  levels 
within  the  beam  region,  and  de-excited  molecules  diffuse  out,  away  from  the  "heap"  of  excess 
de-excited  molecules.  We  treat  all  molecules  in  the  upper  levels  as  a  single  species  (see  Sec.  V 
for  further  details),  diffusing  in  a  background  of  CO^  molecules  in  all  other  levels,  of  N^,  He, 
and  H2;  the  same  is  done  for  the  de-excited  molecules.  The  problem  is  thus  reduced  to  one  of 
diffusion  in  a  6-species  gas:  (a)  the  excited  CO^  molecules,  (b)  the  de-excited  CO^  molecules, 

(c)  CO^  molecules  in  all  other  levels,  (d)  N^,  (e)  He,  and  (f)  H^. 

Consider  a  system  at  constant  temperature  and  constant  total  pressure.  It  may  be  stated 
phenomenologically'1  that  the  density  gradient  of  species  j  is  equal  to  the  weighted  sum  of  the 
diffusion  currents  of  all  species 


v"j  *  2  'A 

k 


(IV-1) 


Onsager's  relations  imply  d.,  =  d,  ..  We  prefer  to  write  the  above  equations  using  different 
symbols  for  the  constants  because  they  are  related  more  directly  to  measured  or  computed  dif¬ 
fusion  constants  in  the  literature.  Also,  when  written  in  the  new  format,  the  equations  automat¬ 
ically  imply  constancy  of  total  pressure  2  Vn.  =  0,  as  appropriate  for  the  present  case.  Then, 

j  1 


Vn 


k^j 


C,  -C. 

Js _ i 

2D.. 

Jk 


D.. 

jk 


D. 


kj 


( IV-2 ) 


The  fact  that  not  all  of  Eqs.  (IV-1)  are  independent  is  displayed  explicitly  in  the  form  of  Eq.  (IV-2). 
Indeed,  if  the  total  pressure  is  to  remain  constant  (note  that  the  temperature  is  constant  and 


*We  gratefully  acknowledge  the  assistance  of  Professor  W.P.  Allis  of  M.I.T.  in  formulating  this  problem. 
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equal  for  all  species),  2  Vn.  =  0,  as  explicitly  shown  in  Eq.  (IV-2)  but  only  implied  by  Eq.  ( IV - 1 ) 

j  3 

through  proper  constraints  imposed  on  the  dj^'s.  Equation  (IV-2)  further  displays  the  physically 
reasonable  fact  that  density  gradients  are  due  to  relative  particle  currents. 

Consider  two  species  as  an  example.  From  Eq.  (IV-2),  it  follows  that 

V(nd  +  n2)  =  0  ( IV —3 ) 

because  =  D^.  This  is  merely  the  condition  that  the  total  pressure  remain  a  constant.  Be¬ 
cause  the  total  particle  current  must  be  zero  to  maintain  the  constant  total  pressure,  =  0. 

It  follows  that 

=  -D12Vni  .  (IV-4) 

This  expression  identifies  as  the  diffusion  coefficient  of  species  1  in  species  2,  for  which 
experimental  data  and  theoretical  expressions  are  available. 

Next,  we  note  that  no  sources  or  sinks  of  the  background  gases  are  present  within  the  dis¬ 
charge.  Therefore,  in  the  steady  state, 

Ck=0  (IV- 5) 

for  k  1,  2,  where  we  identify  with  k  =  1  the  excited  and  with  k  =  2  the  de-excited  CO^  molecules. 

The  current  of  the  excited  molecules  must  balance  the  equal  and  opposite  current  of  de- 
excited  molecules 


Thus, 


We  find  for  the  effective  diffusion  constant  of  the  excited  CO^  molecules 


(IV-6) 


(IV-7) 
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(IV-8) 


Jeans*  * 
molecular  diameters  <j 


cites  a  formula  for  the  diffusion  constants  D^k  in  terms  of  the  molecular  speeds  Cj  and 


3 


Defining 


Sjk  =  2  (ffj  +  V 


CIV-9) 


we  have 


jk 


1 

37^.2 


(IV-10) 


where  v  -  v  .  +  v,  is  the  net  particle  density.  Table  II  gives  the  parameters  used  (taken  from 
3  ^ 

Ref.  12)  and  the  diffusion  constants  D.,  obtained.  The  effective  diffusion  constant  computed  is 
2  -1 

D  ^  46  cm  sec 


11 
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On  the  basis  of  the  foregoing  results,  we  may  estimate  how  important  a  role  diffusion  may 
play  in  determining  gain  saturation.  As  pointed  out  before,  saturation  sets  in  when  the  rate  at 
which  the  optical  field  depletes  the  upper  level  becomes  comparable  to  the  rates  of  relaxation 
of  the  upper  level.  These  relaxation  rates  are  collisional  in  nature,  when  diffusion  is  negligible, 
but  may  become  (in  effect)  enhanced  if  the  particles  diffuse  into,  and  out  of,  the  beam  within  a 
time  comparable  to  the  inverse  collisional  relaxation  rates. 

In  Sec.  Ill,  wc  found  that  the  effective  collisional  relaxation  rate  of  the  upper  level  is  larger 
than  the  actual  rate,  if  the  upper  laser  level  is  coupled  to  other  levels.  The  question  then  is 
what  collisional  relaxation  rate  should  be  employed  to  obtain  an  estimate  of  the  importance  of  dif¬ 
fusion  effects  —  the  actual  or  the  effective  collisional  relaxation  rate?  This  can  be  answered  by 
considering  the  physical  reason  for  the  enhancement  of  the  collisional  over  the  actual  relaxation 
rate  in  a  system  of  M  tightly  coupled  levels,  with  y.U  =  yU  (i  =  1,2,...,  M).  Particles  can  leave 
the  lasing  level  into  the  "general  background"  i.e.,  into  all  the  levels  not  explicitly  included  in 
the  rate  equation  ( 1 1- 1 ),  either  directly  (at  the  rate  yU)  or  by  first  cross -relaxing  into  any  one  of 
the  M-levels,  and  then  leaving  into  the  general  background  at  the  rate  yu.  A  particle  has  effec¬ 
tively  M  ways  of  leaving  the  system  of  upper  levels,  hence  y^  =  My11.  The  saturation  param¬ 
eter  is  raised  as  the  rate  at  which  the  particles  leave  the  upper  level  is  increased.  We  have 
pointed  out  before  that  the  increased  power  density  obtainable  from  such  a  system  (with  given 
small  signal  gain)  is  achieved  because  the  rate  at  which  particles  enter  the  upper  level  is  also 
raised  accordingly. 

Diffusion  is  another  means  by  which  particles  can  enter  (or  leave)  interaction  with  the  laser 
field.  In  a  system  of  M  tightly  coupled  levels,  all  these  levels  diffuse.  Diffusion  effects  will  be 
appreciable  when  the  rate  of  diffusion  of  each  of  these  levels  becomes  comparable  to  the  relaxa¬ 
tion  rate  of  each  of  these  levels,  i.e.,  comparable  to  yU,  not  y^. 

Hence,  when  computing  the  contribution  of  diffusion  to  any  one  particular  group  of  excited  or 
de-excited  CO^  molecules,  we  should  compare  the  rate  at  which  diffusion  supplies  or  extracts 
molecules,  with  the  relaxation  rate  of  a  single  level,  and  not  y^.^. 

We  may  quickly  estimate  the  importance  of  diffusion  effects  by  studying  the  random  walk  of 
an  excited  CC>2  molecule  in  the  background  of  CO^,  N^,  He,  and  H^  molecules  (as  used  exper¬ 
imentally).  By  using  the  effective  diffusion  constant  D  =  46  cm^  sec’\  a  C09  particle  diffuses 
a  distance  d  =  \l 2Dt  =  2.8  mm  within  a  time  t  =  10  3  sec  corresponding  to  the  inverse  relaxation 
rate  l/yu. 

V.  A  THEORY  OF  DIFFUSION  IN  THE  PRESENCE  OF  GAIN  SATURATION 

From  the  foregoing,  it  is  reasonable  to  expect  that  diffusion  plays  an  important  role  in  de¬ 
termining  saturation.  We  now  present  a  simplified  theory  designed  to  evaluate  the  gain  as  a  func¬ 
tion  of  beam  geometry  and  input  intensity,  and  hence  the  effective  saturation  parameter  I  as 
influenced  by  diffusion. 

To  include  diffusion  effects  in  the  rate  equations  of  a  system,  a  particle  current  can  be 
defined  which  is  proportional  to  the  gradient  of  the  number  density,  with  the  proportionality  fac¬ 
tor  being  the  diffusion  constant  D 

C  =  -DVN  .  (V-l) 
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Then,  from  the  conservation  of  mass. 


v-  cp 
dv2n 


_f9Nl 

1  at  J  diffusion 

f9N, 

1  8t  J  diffusion 


Formally,  let 


-DV 


2 


and 


( V-2) 
(V-3) 


ft  f.  u  u 

Then,  Eqs.  ( II  —  4 )  still  hold  formally,  with  y .  replaced  by  I\  ,  and  y.  by  I\  .  Actually,  the  prob¬ 
lem  now  involves  a  set  of  coupled  partial  differential  equations.  We  can  reduce  this  set  to  two 
coupled  equations  in  the  limit  of  fast  rotational  relaxation  by  the  following  heuristic  argument: 

In  the  absence  of  diffusion,  for  tight  coupling,  the  saturation  intensity  is  that  of  a  two-level  sys¬ 
tem,  with  y11  represented  by  2  y  U(N  e/N^ ),  and  y  represented  by  £  y/(n.e/n^).  If  diffusion  is 

i  u  i  i  i 

treated  formally  by  replacing  y^  by  I\  ,  and  y.  by  I\  ,  the  expression  for  the  saturation  inten¬ 
sity  is  that  which  would  have  been  obtained  from  a  two-level  system  with  the  diffusion  constant  D 
replaced  by 


N. 

D'  =  D  — 

U  Nf 

1 


for  the  upper  level,  and 


Di 


D 


n. 

i 


for  the  lower  level.  We  may  use  this  argument  to  replace  the  set  of  coupled  partial  differential 
equations  by  two  such  equations.  Further,  for  simplicity,  we  use  the  same  diffusion  constant 
for  both  levels 


M 

D'-D  2 

i=  1 


Further,  using  the  symbols 
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we  obtain 


/  g  V 

D'V2Nd  -ruN1  -  w(n4  -  ^ni)=  -R*  (V -4) 

and 

D'V2^  -y^  +  w(n4  nj  =  -R*  .  (V-5) 

Solution  of  Eqs.  (V-4)  and  (V-5)  for  a  Gaussian  beam  is  quite  difficult.  However,  a  relatively 
straightforward  solution  can  be  obtained  if  the  beam  has  no  radial  variation  in  intensity.  Allow¬ 
ing  W  to  be  constant  across  the  beam  and  setting  it  zero  outside,  we  obtain  by  combining 
Eqs.  (V-4)  and  (V-5), 

(D’V2  _yu  -W)  (D'V2  -yt)  -  W(D’V2  -yu)J  =  TuR/+  W(r“  +  R*)  .  (V-6) 


Assuming  a  modified  Bessel  function  for  the  homogeneous  solution  leads  to 


Y  UR/  +W(RdU  +  R/) 

(gu  \ 

TuY,  +W(-  Yu  +  Yi) 


+  ^+I0(<P+r)  +  ’PJ0^.r'> 


(V-7) 


where 


and  ip  and  ip  are  constants  to  be  determined  from  the  boundary  conditions.  Since 


Ni  =  -w  [(D'v2  -  |j  w)  ni  +  R/ 


the  population  inversion  of  the  medium  inside  the  beam  can  be  written  as 
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Outside  the  beam,  W  becomes  zero  and 
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where  and  are  constants  to  be  obtained  from  the  boundary  conditions.  We  require  conti¬ 
nuity  of  population  density  and  particle  current  at  the  beam  boundary  r  .  Defining 
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the  boundary  conditions  can  be  written  in  matrix  form  as 
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It  is  also  convenient  to  define 


G(a,b)  = 


I  (<p  r  ) 
oXYa  o' 


(p  I  Aw  r  )  + 

a  l'v  a  o 


~Ko(J^ro) 
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( V — 1 7 ) 


where  a  assumes  the  "values"  +  and  — ,  and  b  assumes  the  "values"  t  and  u.  Then, 

(ir--^u)Ki  (n/¥ ro) +  a-(^7 Ki(J5^)J¥e<-’u> 


>P  = 
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a+g(+,  u)  e(-  i )  -  A  _e(-  u)  e(+,  i ) 


(V-18) 


(V-19) 


Equation  (V-10)  shows  that  diffusion  causes  both  the  small -signal  gain  and  the  saturation 
parameter  to  become  functions  of  position  in  the  beam.  However,  since  most  gain  expressions 
in  the  literature  assume  no  radial  variation  in  saturation  parameter,  some  sort  of  "average M 
value  is  desirable  for  comparison  of  results.  Also,  a  Gaussian  beam  would  probably  have  a 
much  different  radial  variation  than  the  uniform  beam,  again  implying  that  some  sort  of  averag¬ 
ing  should  be  done.  With  this  in  mind,  we  can  define  an  average  gain  a 
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Using  Eqs.  (V-9),  (V-18),  and  (V-19)  in  (V-20)  gives 
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As  can  be  seen,  a  closed-form  expression  for  the  saturation  parameter  I  cannot  be  defined  in 

the  usual  way  from  Eq.  (V-21).  However,  I  can  be  found  graphically.  Taking  I  as  the  intensity 

s  s 

at  one-half  the  small-signal  gain  gives  the  saturation  parameter  for  a  two-level  system.  The 
multilevel  saturation  parameters  are  plotted  as  a  function  of  beam  radius  in  Fig.  1  which  is  a 
plot  of  I  obtained  from  Eq.(V-21)  as  that  intensity  for  which  the  gain  drops  to  half  its  small- 
signal  value.  A  diffusion  constant  of  46  cm  /see  as  computed  in  Sec.  IV  for  one  diffusion  species 


Fig.  1.  Computed  variation  af  saturation 
parameter  with  and  without  diffusion. 


0  2  4  6 

rQ  (mm) 


has  been  used.  The  number  of  equivalent  levels  participating  affects  the  evaluation  of  the  sat¬ 
uration  parameter  only  in  one  way:  when  we  search  for  the  W(Ig),  which  reduces  the  gain  to 
half  its  value  at  zero  intensity,  we  find  that  W(Ig)  is  proportional  to  the  relaxation  rate  (if 
both  yu  and  y ^  are  increased  by  the  same  factor).  But  this  implies  that  I  is  proportional  to  the 
relaxation  rate  and  hence  to  the  number  of  equivalent  levels.  Figure  1  has  been  plotted  for  14 
levels  as  computed  in  See.  IV.  We  have  used  the  parameters  given  in  Table  III  below. 


TABLE  III 

Parameter 

Val  ue 

Reference 

t 

SP 

5.2  sec 

13,  14 

1 

r  u 

u 

10  ^ 5ec 

15 

1 

rT< 

»T 

U 

- 

Ar 

43  MHz 

16 

18 


VI.  EXPERIMENT 


Before  discussing  the  details  of  the  experiment,  we  want  to  develop  the  theory  used  to  ex¬ 
tract  from  the  measurements  the  value  of  the  saturation  parameter  as  a  function  of  geometry, 
a  parameter  that  should  be  independent  of  geometry  if  diffusion  effects  were  negligible.  In  or¬ 
der  to  avoid  the  very  difficult  analysis  entailed  by  inclusion  of  diffusion,  we  treat  the  system  as 
if  it  were  unaffected  by  diffusion,  predict  the  curve  of  gain  vs  intensity  on  this  basis,  and  then 
match  the  experimental  results  to  this  curve.  For  different  experimental  beam  diameters,  we 
then  find  different  saturation  parameters.  We  take  this  as  proof  of  the  importance  of  diffusion 
in  determining  saturation. 

Consider  first  the  expression  for  the  gain  constant  a  in  the  limit  of  small  intensity  I.  Re¬ 
gardless  of  the  model  of  the  system  (whether  homogeneously  or  inhomogeneously  broadened), 
as  long  as  the  analysis  is  based  on  a  rate  equation  without  inclusion  of  diffusion  effects,  the  de¬ 
pendence  is 


n  -  n 

o 


(VI-1) 


where  the  above  equation  can  serve  as  a  definition  of  the  parameter  I  [compare  Eq.  (11-12)). 

Now  assume  that  diffraction  is  unimportant  so  that  ray  optics  can  be  applied.  Further  as¬ 
sume  that  the  optical  beam  rays  are  parallel,  a  situation  which  is  reasonably  well  approximated 
in  the  experiment  (see  Table  IV).  Then,  at  any  radial  distance  r,  I(r,  z)  is  given  by 


1 

I  dz 


a 


(VI-2) 


TABLE  IV 


r 

max 

(mm) 

r 

(mm) 

r  . 
mm 

(mm) 

di 

(m) 

d2 

(m) 

R 

c 

(m) 

r 

meas 

(mm) 

r 

comp 

(mm) 

1.1 

0.9 

0.7 

4.73 

2.19 

3.12 

2.6 

2.5 

1.2 

1.1 

1.1 

3.34 

2.36 

3.12 

2.0 

1.9 

1.8 

1.8 

1.8 

3.80 

4.16 

4.66 

2.0 

2.0 

2.5 

2.5 

2.5 

2.46 

2.62 

4.66 

2.6 

2.6 

r  =  maximum  l/e  radius  of  intensity 
max 


r  =  average  l/e  radius  of  intensity 

r  .  =  minimum  l/e  radius  of  intensity 
min  7  7 

d^  =  oscillator  output  mirror  to  spherical  mirror  separation 

d 2  =  amplifier  center  to  spherical  mirror  separation 

=  spherical  mirror  radius  of  curvature 

r  =  measured  l/e  radius  of  intensity  1  m  from  amplifier  center 
meas  7 

fcom  =  comPu*ec^  l/e  ra^'us  intensity  1  m  from  amplifier  center 
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Integrating  Eq.  (VI-2)  results  in 


G  I  I. 
o  o  in 


'out  -  I0  +  (G0-l)Iin 


a  L 


where  G  =  e  °  ,  and  L  is  the  length  of  the  amplifier.  Experimentally,  we  measure 


out 


■1 


I  ,27rrdr 
out 


(VI-3) 


(VI-4) 


and 


P. 

in 


I  hr 


27rrdr 


(VI-5) 


Integration  of  the  above  for  an  assumed  Gaussian  intensity  distribution 

P.  ,  2,-2. 

I.  (r)=  -!£  e-(r  /r  > 
in  -2 

7rr 


(VI-6) 


gives 


G  = 


P  , 
out  _  - 

“p-  -  Go 


1  - 


G  -  1  P. 
o  in 


~  2  - 
7rr  I 

o 


(VI-7) 


Therefore,  when  plotting  gain  G  as  a  function  of  P.  /vv  ,  we  would  obtain  the  same  straight  line 
~2  m 

for  all  values  of  r  if  I  were  independent  of  the  beam  diameter.  Such  a  plot  then  serves  as  a 
means  of  determining  I  .  Suppose  next  that  a  is  not  given  simply  as  Eq.  (VI-1),  but  has  the  gen¬ 
eral  form 


a  =  f(a  ,  I,  I  ) 
o  o 


(VI-8) 


where  f  is  some  decreasing  function  of  I  whose  initial  value  is  n  and  whose  initial  slope  against 
I  is  a  /i  .  Integrating  Eq.(VI-8)  results  in  a  new  function 


I 


out 


g( O'  ,  I .  ,  I  ,  L) 
o  in  o 


(VI-9) 


where  L  is  the  length  of  the  amplifier  as  before.  Evaluating  the  gain  then  yields 

2  /-2X 
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(VI-10) 


Thus,  we  obtain  for  gain  G  a  function  of  I  ,  (P.n/7rr2),  aQt  and  L.  A  plot  of  gain  vs  P.n/irr 
for  varying  average  beam  radii  r  should  give  a  single  plot,  if  no  other  parameters  are  varied 
in  the  amplifier.  Deviations  from  such  a  behavior  are  an  indication  that  beam  geometry  plays 
a  role. 
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Fig.  2.  Computed  gain  variation  without  diffusion  for  different  values 
of  saturation  parameter. 

—2 

To  indicate  the  dependence  of  gain  upon  the  parameter  (pin/^r  ),  Fig.  2  shows  plots  of  G 
for  varying  Ig  =  I  assuming  homogeneous  broadening.  The  straight-line  slopes  correspond  to 
Eq.  (VI-7). 

Now  let  us  turn  to  the  experiments  designed  to  measure  the  influence  of  diffusion  upon  sat¬ 
uration.  Figure  3  shows  the  equipment  and  experimental  layout.  The  signal  source  was  a  stable, 
sealed-off  oscillator  operating  in  a  single  TEM^^  mode  with  up  to  15-W  output  power.  The 
cavity  configuration  of  the  oscillator  was  semi-confocal  with  mirrors  spaced  451.3  cm  apart,  and 
the  output  beam  had  a  l/e  radius  of  4.6  mm.  In  this  report,  beam  radius  is  defined  by  the  radial 
distance  at  which  the  intensity  of  the  beam  (not  the  electric  field)  falls  to  l/e  times  its  value  at 
the  beam  center.  The  design  and  stability  of  the  laser  oscillator  has  been  discussed  elsewhere.*^ 


Fig.  3.  Experimental  setup. 
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In  order  to  adjust  the  average  beam  radius  in  the  amplifier  and  keep  the  beam  within  the 
perimeter  of  the  3  X  8-ft  optical  table,  the  output  of  the  laser  was  usually  made  to  follow  a  path 
between  four  mirrors,  one  of  which  was  curved.  In  two  cases,  the  short  path  lengths  necessary 
for  matching  allowed  the  use  of  only  three  mirrors.  The  beam  was  directed  into  the  amplifier 
in  a  confocal  configuration  with  the  beam  waist  at  the  center.  The  average  beam  diameter  in  the 
amplifier  was  controlled  by  adjusting  the  path  lengths  between  the  oscillator,  the  spherical  mir¬ 
ror,  and  the  amplifier,  and  by  choosing  the  proper  radius  of  curvature  from  available  matching 

mirrors.  Path  lengths  and  mirror  curvatures  for  the  various  average  beam  radii  are  given  in 

1 8 

Table  II.  Path  lengths  were  calculated  with  the  mode-matching  formulas  of  Kogelnik  and  Li. 
Beam  radii  were  checked  by  centering  a  4-mm-diameter  iris  in  the  beam  100  cm  from  the  am¬ 
plifier  center  and  using  the  relation 


where 

r^  =  l/e  radius  of  the  intensity 
r1  radius  of  iris 
P1  =  power  with  iris  in 
P  =  power  with  iris  out. 


(VI-11) 


Equation  (VI-11)  was  obtained  by  integrating  the  beam  intensity  over  the  area  of  the  iris  to  find 
the  power  transmitted  through  the  iris.  The  measured  beam  radii  were  compared  with  computed 
values  by  using  the  Gaussian  beam  formulas  and  assuming  correct  waist  diameter  and  location. 
Maximum  deviation  of  the  measured  radii  was  about  7  percent,  and  was  probably  due  to  astig¬ 
matism  caused  by  the  use  of  a  spherical  mirror  in  an  off-normal  incidence. 

The  water-cooled  amplifier  had  an  inside  diameter  of  1.8  cm,  and  70  cm  of  active  length.  A 
sealed-off  amplifier  was  used  because  pumping  speeds  are,  in  general,  comparable  to  diffusion 
speeds  and  would  tend  to  contribute  in  an  incontrollable  way  to  the  diffusion  effects.  A  relatively 
large  inside  diameter  (4  to  10  times  larger  than  the  beam  diameters)  was  chosen  to  minimize 
wall  effects.  The  amplifier  was  filled  with  a  mixture  of  7  torrs  He,  1.4  torrs  CO^,  1.7  torrs 
N^,  and  0.2  torr  II^.  Discharge  was  from  a  center  cathode  to  anodes  at  each  end.  Optimum 
cathode  current  was  found  to  be  26  mA.  The  beam  entered  the  amplifier  through  NaCl  windows 
tipped  7°  to  the  beam  axis  to  reduce  the  possibility  of  feedback  effects.  Insertion  loss  of  the 
amplifier  with  the  discharge  off  was  about  18  percent.  Nearly  all  this  loss  was  due  to  reflection, 
scatter,  and  absorption  by  the  two  windows. 

The  total  power  output  of  the  amplifier  output  was  measured  without  an  iris  by  using  a  Co¬ 
herent  Radiation  Laboratory  thermopile  capable  of  directly  measuring  the  output  power.  Input 
power  to  the  amplifier  was  varied  by  a  combination  of  controlling  the  oscillator  current,  using 
only  one  anode  of  the  oscillator,  and  replacing  totally  reflecting  gold-coated  mirrors  with  par¬ 
tially  transmitting  dielectric  mirrors. 

Since  there  was  no  significant  difference  in  insertion  loss  with  the  amplifier  completely  evac¬ 
uated  and  with  the  amplifier  filled  with  the  gas  mixture,  absorption  in  the  CO^  was  neglected  in 
all  measurements.  Power  gain  was  taken  as  the  ratio  of  output  power  with  the  discharge  on,  to 
output  power  with  the  discharge  off.  The  sum  of  the  output  power  with  the  discharge  off  and  the 
power  lost  in  the  exit  window  was  used  as  the  input  power  to  the  amplifying  medium. 
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Fig.  4.  Experimental  resul ts  af  gain  variation  far  various  average  beam 
radii  as  a  function  af  input  power  density  in  watts  per  square  centimeter. 

Gain  measurements  were  made  for  four  different  average  beam  radii  between  0.9  and  2.5mm 
(measured  to  the  l/e  point  of  the  intensity).  Figure  4  shows  power  gain  as  a  function  of  input  in¬ 
tensity  for  the  various  beam  radii  used.  The  saturation  parameters  were  obtained  from  the  plots 
of  Fig.  4,  using  the  initial  slopes  as  indicated  by  Eq.  (VI-2).  Note  the  constancy  of  the  small- 
signal  gain  which  indicates  the  relative  unimportance  of  wall  effects. 

Figure  5  shows  the  saturation  parameter  as  a 
function  of  beam  radius  and  is  a  summary  of  the 
experimental  results.  The  measured  saturation 
parameter  varied  from  25  to  97  W/cm  ,  and  was  a 
monotonically  decreasing  function  of  beam  radius. 

In  conclusion,  experimental  evidence  indicates 
that  diffusion  effects  can  play  an  important  role  in 
determining  the  saturation  parameter  in  CO^  la¬ 
sers.  The  experimentally  measured  variation  of 
saturation  parameter  agrees  qualitatively  with  a 
curve  obtained  from  a  simple  theoretical  treatment 
of  diffusion  effects.  Not  too  much  importance  can 
be  given  to  the  quantitative  results  of  the  diffusion 
theory  which  has  been  worked  out  for  an  optical 
beam  of  a  "square”  rather  than  a  Gaussian  cross 
section.  If  we  attempted  to  match  the  theory  to  the 
experiment,  a  value  of  50  (rather  than  14)  equiv¬ 
alent  levels  would  give  better  agreement  with  ex¬ 


Fig.  5.  Experimentally  measured  saturation 
parameter  as  a  function  af  average  beam 
radius. 


periments.  This  discrepancy  can  be  due,  in  part,  to  the  difference  in  the  geometry  analyzed  from 
that  used  in  the  experiment;  it  could  also  be  due  to  an  underestimate  of  the  equivalent  number  of 
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levels  when  the  number  14  was  found  in  Sec.  IV.  The  vibrational  temperature  may  be  consider¬ 
ably  higher  than  the  rotational  temperature,  resulting  in  a  larger  number  of  effective  levels. 
Further,  in  addition  to  the  CO^  levels,  some  levels  are  also  tightly  coupled  into  the  lasing 
level,  and  this  coupling  may  also  raise  the  number  of  effective  levels.  A  more  rigorous  treat¬ 
ment  of  the  problem  is  yet  to  be  done  and  will  lead,  undoubtedly,  to  more  refined  experiments. 
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are  comparatively  low  so  that  one  would  expect  spatial  diffusion  to  play  an  important  role  in  determin¬ 
ing  saturation.  A  theory  of  diffusion  is  carried  through  and  it  is  shown  that  diffusion  effects  can  indeed 
be  important  for  optical  beam  diameters  of  a  few  millimeters.  Finally,  experiments  on  a  sealed-off  CO2 
laser  oscillator  and  amplifier  system  are  reported.  Amplifier  gain  is  measured  for  four  different 
beam  radii.  The  "equivalent”  saturation  parameter  derived  from  the  measurements  decreased  monoton  - 
ically  from  97  to  25  W/cm^  as  the  average  input  beam  radius  increased  from  0.  9  to  2.  5  mm  in  the  9-mm 
radius  discharge  tube  of  the  amplifier. 
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